Signatures of the Origin of High-Energy Cosmic Rays in Cosmological 

Gamma-Ray Bursts 



Jordi Miralda-Escude and Eli Waxman 
Institute for Advanced Study, Princeton, NJ 08540 



> 

(N 



o 

:^ 

o 

u 



ABSTRACT 

We derive observational consequences of the hypothesis that cosmic rays (CR's) 
of energy > lO^^eV originate in the same cosmological objects producing gamma-ray 
bursts (GRB's). Inter-galactic magnetic fields ^ 10~^^G are required in this model 
to allow CR's to be observed continuously in time by producing energy dependent 
delays in the CR arrival times. This results in individual CR sources having very 
narrow observed spectra, since at any given time only those CR's having a fixed time 
delay are observed. Thus, the brightest CR sources should be different at different 
energies. The average number of sources contributing to the total CR flux decreases 
with energy much more rapidly than in a model of steady CR sources, dropping to 
one at Scrit — 2 x 10^'^ eV with very weak sensitivity to the inter-galactic magnetic 
field strength. Below -Ecrit, a very large number of sources is expected, consistent 
with observations. Above i^crit, a source may be observed with a flux considerably 
higher than the time-averaged CR flux from all sources, if a nearby GRB occurred 
recently. If such a source is present, its narrow spectrum may produce a "gap" in the 
total spectrum. These signatures should be detectable by the planned "Auger" CR 
experiment. 
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Introduction 



The sources of gamma ray bursts (GRB's) and of cosmic rays (CR's) with energy E > lO^^eV 
are unknown. In particular, most of the sources of cosmic rays that have been proposed have 
difficulties in accelerating CR's up to the highest observed energies ( piUas 1984 , Cronin 1992| ). 
Recent gamma ray and cosmic ray observations give increasing evidence that both phenomena are 
of cosmological origin [see [Meegan et al. 1992| , paczyhski 19921 [Piran 1992| for GRB's; [Bird et al 



1994| , [Yoshida et al. 1995| , [Waxman 1995b| for CR's]. Although the source of GRB's is unknown, 
their observational characteristics impose strong constraints on the physical conditions in the 
7-ray emitting region ( piran 1994 , Meszaros 1995| ), which imply that protons may be accelerated 
in the 7-ray emitting region to energies 

1920 _ io2ieV ( IWaxman 1995a| [Vietri 1995| ). In addition, 
the average rate (over volume and time) at which energy is emitted as 7-rays by GRB's and in 
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CR's above lO^^eV in the cosmological scenario is, remarkably, comparable ( [Waxman 1995a| ,b). 
These two facts suggest that GRB's and high-energy CR's may have a common origin. 

An essential ingredient of a bursting model for cosmic rays is the time delay due to inter- 
galactic magnetic fields. The energy of the most energetic CR detected by the Fly's Eye experiment 
IS m excess of 2 X lO^OeV ( pird et al. 199^ , and that of the most energetic AGASA event is above 
lO^'^eV ( Payashida et al. 1994| ). On a cosmological scale, the distance traveled by such energetic 
particles is small: < lOOMpc for the AGASA event, < 50Mpc for the Fly's Eye event ( [Aharonian 



fc Cronin 1994 ). Thus, the detection of these events over a ~ 5yr period can be reconciled with 
the rate of nearby GRB's (~ 1 per 50yr in the field of view of the CR experiments out to lOOMpc 
in a standard cosmological scenario; e.g., [Cohen &: Piran 1994 ) only if there is a large dispersion 



in the arrival time of protons produced in a single burst. The required dispersion, > 50yr for 
10^*^ eV proton, may be produced by inter-galactic magnetic fields (Waxman 1995a, [Waxman & 



Coppi 1995[ ). The deflection angle for a proton propagating a distance D in a magnetic field B 



with coherence length A is ~ 0.05° (D/A)^/^ (A/lOMpc) (5/10^^^ G) (-E/1020eV)"^ and the 
induced time delay is t{E) ~ 103yr(Z)/100Mpc)2 (A/lOMpc) (S/IO"" G)^ {E / eVy^ . Since 
the time delay is energy dependent, the large spread in proton energy, induced by random energy 
loss, results in a time broadening of the CR pulse over a time ~ t{E). Thus, the CR's from a 
single burst can be received on Earth over a long time interval. Nevertheless, since the angular 
deflection is small, the individual sources are still detectable by measuring the arrival directions. 

In this Letter, we examine the characteristics of CR sources that should be expected in a 
bursting source model, as a function of the time delay, dependent on the magnetic field. We find 
that there are characteristic signatures for such a model, which would allow to distinguish it from 
a scenario where the CR sources are steady, i.e., where the sources emit a constant flux on a time 
scale longer than the time delay of the lowest energy CR's that are relevant. In §2.1 we give a 
qualitative description of the bursting model properties, using an approximate analytic approach. 
In §2.2 we present Monte-Carlo simulations that demonstrate the properties discussed in §2.1. The 
various tests of the bursting model and the implications for future high energy CR experiments 
are discussed in §3. 



2. Characteristics of Cosmic Ray Bursts Sources 

We consider a Cosmic Ray Burst (CRB) taking place at a distance D from us, where a total 
number of protons np{E)dE of energy E is emitted at a single instant in time. The CR's arrive 
with a time delay t, relative to gamma-rays, with a probability density p[t/T{D,E)\(l[t/T{D,Ey\, 
where t{D,E) oc D'^/E'^ is the characteristic time delay. The time delay t at a fixed energy and 
distance varies randomly due to two effects. First, the magnetic field along a trajectory should 
have random variations; for example, if the magnetic fields originate in galaxies and are later 
ejected to the inter-galactic medium, the field strength along a trajectory should vary depending 
on the impact parameter to individual galaxies. In the absence of energy losses, the bursting 
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source would produce a number of cosmic ray "images" , and the cosmic rays in each image would 
be of a single energy which would decrease with time as However, the random nature of 

the energy loss of a cosmic ray eliminates these images, and simply introduces a dispersion in the 
arrival times and arrival directions at a fixed energy. In general, the dispersion in arrival times t 
will be of order t{D,E) ( |Waxman fc Coppi 1995| ). 



2.1. Analytic Model 

We now perform a simple analytic calculation of the number of CR sources that should be 
seen at each energy and flux in a CRB model. For this purpose, we approximate the effect of 
energy losses as being negligible when CR's come from a distance D < Dc{E), and eliminating 
all cosmic rays coming from D > Dc{E) (for E < 10'^'' eV, this approximation is quite good and 
Dc{E) corresponds to the distance where the initial proton energy necessary to have an observed 
energy E, after losses to electron-positron production, exceeds the threshold for pion production). 
We also assume that the sources are observed only during a time t{D,E) with a constant flux 

""P^^^ _ n,{E)D,{Ef 
^ A7rD'^T{D,E) ATiD^TciE) ' 

where Tc{E) = t{Dc{E), E). If the rate per unit volume of CRB's is all emitting the 
same np{E), then the average number of bursts at distance D observed at any time is 
n{D,E) dD = Atti'Tc{E)[D'^ / Dc{E)'^] dD, giving a number of bursts at a given observed flux 



n{F,E)dF = ^vDc{EfTc{E) ' 



F 



F 



(2) 



where Fc{E) = np{E) /[ATTDc{EfTc{E)]. The flux Fc{E) is the minimum flux observed for the 
sources. In our simplified model, the number of sources drops to zero abruptly at Fc{E) owing 
to the assumed distance cutoff Dc{E) and the "top-hat" time profile. In reality, there should be 
a smooth turnover at Fc{E) of the number of CRB sources from the —5/4 power-law slope at 
the bright end. This result for bursting sources is in contrast to the usual —3/2 Euclidean slope, 
which applies for steady sources of cosmic rays. 



The total average number of sources above flux F is 



Attv 



N{F,E) = —D,{EYt,{E) 



Fc{E) 



5/4 



(3) 



F 

and the average background flux resulting from all the sources is 

B{E) = AT:uD,{EfT,{E)F,{E) = vnp{E)D,{E) . (4) 

The background flux is dominated by sources with flux near Fc{E), although the contribution 
from brighter sources decreases only as 
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As the cosmic ray energy is increased, the average number of bursts observed above the 
turnover flux Fc{E) decreases, and there is a critical energy -Ecrit where this average number of 
sources equals unity: 



(5) 



We can write the average number of sources in terms of -Ecrit as 



N,{E) = N[F,{E),E] 



E, 



crit 



Dc{E) 

Dc{Ecr\t) 



(6) 



The number of sources Nc drops rapidly with energy, due to the strong dependence 
on the decreasing cutoff distance Dc{E). The drop is especially rapid near lO^'^eV, where 
Dc{E) decreases quickly (see Fig. 2 below). Therefore, for E < ii^crit, the number of 
sources contributing to the flux is very large, and the total number of CR's received at any 
given time is near the average background B{E). The brightest source has a typical flux 
Fi{E) ~ F,{E)N,{E)^/^ = [B{E)/5]iE/E„itf/^[D,{E)/D,iE,,,t)]-\ although there is a 
probability to observe a source with F > Fi(S), P ~ [F/ Fi{E)]-^/^ . At E > Ecrit, the total 
energy received in CR's will generally be much lower than the average B{E), because there will 
be no burst within a distance Dc{E) having taken place sufficiently recently. The few CR's may 
be the lucky survivors from sources further than Dc{E), or they may have anomalously long 
time-delays as a result of crossing a region of high magnetic field (probably near a galaxy) . There 
is, however, a probability P ~ Nc{E) of seeing one CR source with E > E^rit having a flux 
~ B{E)/Nc{E), or an even brighter one with probability decreasing as F~^^^. 

If the CR sources are steady, then the number of sources decreases with energy only as 
Dc{EY , i.e., much more slowly than predicted by eq. (^). This implies that for a given critical 
energy, the number of bright sources aX, E < -Ecrit predicted by a model of steady sources is much 
larger than that predicted for bursting sources. 

Bursting CR sources should have narrowly peaked energy spectra, and therefore the brightest 
sources should be different at different energies. For example, if a bright source is observed at 
E > -Ecrit) the burst must have taken place recently in order that the high energy cosmic rays are 
just arriving on Earth, so lower energy cosmic rays will not have arrived yet. Typically, there 
will be other brighter sources at E < Ecrit, corresponding to bursts that took place a longer time 
ago, and probably closer to us (since there is a longer time interval available). This is in marked 
contrast to a model of steady state sources, where the brightest source at high energies should also 
be the brightest one at low energies, its fractional contribution to the background decreasing to 
low energy only as Dc{E)^^. 
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2.2. Numerical Results 

We now present the results of a Monte-Carlo simulation of the total number of cosmic rays 
received from CRB's at some fixed time. For each realization we randomly draw the positions 
(distances from Earth) and times at which cosmological CRB's occured, assuming that the CRB's 
are homogeneously distributed standard candles with an average rate v = 2.3 x 10^'^/i'^Mpc~'^yr~^ 
(with h = 0.75) similar to that fitted to the observed flux distribution of CRB's assuming a 



no-evolution standard candle model ( Cohen fc Piran 1994 ). We assume an intrinsic cosmic ray 



generation spectrum np{E) cx E~'^dE, which produces a flux above 2 x lO^^eV consistent with 
the Fly's Eye and AC AS A data (Waxman 1995b). We calculate the change of the spectrum due 
to interaction with the CMB photons in a method similar to that described in Waxman (1995b), 
except that for distances < ISOMpc we do not use the continuous energy loss approximation but 
rather an exact calculation of the energy loss, which includes fluctuations. 

For the probability distribution of the time-delay for a cosmic ray of fixed energy from a source 
at a given distance, we use the form p(r) = Po{t/to)~'^~^ for > tq, p{T)dT = {4po/3){t/tq — 1/4) 
for To > T > ro/4, and p(r) = for r < ro/4, where tq is the characteristic time-delay and is 
proportional to D'^/E'^, and po = [ro(3/8 -|- l/a)]~^. For r < tq the form of the probability 
function approximately matches results obtained by Waxman & Coppi (1995). Cosmic rays 
with large r are the ones that have crossed regions of high magnetic field, and we assume 
there is a power-law distribution of magnetic field strengths, giving a power-law distribution 
of time delays. If the typical field in the inter-galactic medium is B and has coherence length 
A, the typical deflection angle is 6 (x B (DX)^^^ . When intercepting a region with magnetic 
field B' ^ B and coherence length A', the deflection angle is 9' oc B'X' , yielding a time delay 
r' oc t{B'X' /BX)^{X/D). If n is the number density of such regions, the interception probability 
is nvrA'^D, so the index a is q = —log{iTnX''^D)/log[{B''^X''^)/{B'^XD)]. Here, we shall use as 
an example ti = t{D = 80Mpc, E = lO'^'^eV) = 10'^ years, corresponding to A ~ lOMpc and 
B We also take B'/B = 10^ A'/A = 10-^ and n = IQ-^Mpc"^, giving typical 

parameters for spiral galaxies (this leads to q ~ 1, with a weak dependence on distance). 

We have examined a total of 50 realizations. About 70% of these are similar to the realization 
presented in Fig. la: there is no source sufficiently nearby having occurred sufficiently recently, 
so the flux at high energies is below the average. In the other 30%, there are typically one or two 
bright sources dominating at -E > lO^'^eV, as in the realization presented in Fig. lb. A source 
similar to the brightest one in Fig. lb appears only 4% of the time (in this example, the source 
is at z = 0.0056 and occurred 51 years ago); the second brightest source at -E ~ lO^'' eV in Fig. 
lb is more common. The analytic expression (^) for N{F, E) provides a good approximation to 
the numerical results for the number of sources with flux F and spectral peak at E, except that 
sources at high energy are also present at fluxes below Fc{E), coming from CRB's at distances 
higher than Dc{E) for which some high energy cosmic rays still survive. The spectral shape of the 
individual sources is determined by the time-delay probability distribution we have assumed, and 
is slightly modified by the interaction with the microwave background (this is the reason why the 
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shape of the spectra of different sources varies). 

Fig. 2 shows Nc{E), calculated from the average background and equations (^) and 
(^. Since our numerical model does not assume a sharp cutoff in the flux distribution, as 
we did above in the analytical model, Nc{E) is here an indication of the number of sources 
above the turnover flux, which dominate the contribution to the average background. All the 
characteristics of the sources depend on the CRB rate v and on the characteristic time delay 
Ti = t{D = 80Mpc, = 10^'^eV) only through their product vti, or, equivalently, through 
the critical energy i?crit(i^''"i)- For the parameters we have chosen, z/ti = 

lO^^Mpc-^, Fig. 2 

shows that Scrit — 1-4 x 10^'' eV. The dependence of -Ecrit on uti is easily determined from 
Fig. 2, since Nc oc uti (see eqs. P-^) and therefore the curve in Fig. 2 shifts vertically as 
uTi. If > lO^^eV CR's are indeed produced by GRB's, then u is determined by the GRB flux 
distribution. The time delay, however, depends on the unknown properties of the inter-galactic 
magnetic field, ri oc B^X. As mentioned in §1, current data requires ri ^ 50yr, or, equivalently, 
£^crit ^ 10^° eV, which corresponds to BX^/^ ^ lO^^^G Mpc^/^. The current upper limit for 
the inter-galactic magnetic field, BX^/^ < lO'^G Mpc^/^ ( [Kronberg 199^ , [Vallee 1990] ), allows 



a much larger delay, ri < lO^yr. However, the rapid decrease of Nc{E) with energy near 
lO^'^eV, implies that -Ecrit is not very sensitive to vti. Thus, for the range allowed for the 
GRB model, 5 x IQ-'^Mpc^^ < uti < IQ-^Mpc"^, the critical energy is limited to the range 
lO^OeV < Ecrit < 3 X lO^OeV. 



3. Discussion 

We have analyzed a model where > lO^^eV CR's are produced by cosmological sources 
bursting at a rate comparable to GRB's. We have found that, in this model, the average 
number of CR sources contributing to the flux decreases with energy much more rapidly than 
in the case where the CR sources are steady. We have shown that a critical energy exists, 
lO^OeV < Ecrit < 3 X lO^OeV, above which a few sources produce most of the CR's, and that the 
observed spectra of these sources is very narrow: the bright sources at high energy should be 
totally absent in cosmic rays of substantially lower energy, since particles take longer to arrive the 
lower their energy. In contrast, a model of steady sources predicts that the brightest sources at 
high energies should also be the brightest ones at low energies. 

Above i^criti there is a significant probability to observe one source with a flux considerably 
higher than average. If such a source is present, its narrow spectrum may produce a "gap" in the 
overall spectrum, as in Fig. lb. Recently, ^igl et al. (1995) argued that the observation of such an 



energy gap would imply that the sources of > 10^'^ eV CR's are different from the sources at lower 
energy, hinting that the highest energy CR's are produced by the decay of a new type of massive 
particles. We see here that this is not the case when bursting sources are allowed, owing to the 
time variability. If such an energy gap is present, our model predicts that most of the cosmic rays 
above the gap should normally come from one source. 
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If our model is correct, then the Fly's Eye event above 2 x 10^'^eV suggests that we live 
at one of the times when a bright source is present at high energies. However, the absence of 
such a source can not be ruled out, since, for example, the probability to have detected the Fly's 
Eye event in the realization of Fig. la, where no bright source exists, is ~ 3%. The highest 
energy AGASA event might more easily be produced by a common, faint source (like in Fig. la). 
Furthermore, notice that, given that Fly's Eye has detected only one cosmic ray with E > 10^'^ 
eV, we already know that the AGASA cosmic ray had a low probability of being detected; within 
the measurement error, its energy might be not much above lO^'' eV. 

Given the present scarcity of ultra-high energy CR's, no solid conclusions can be drawn. 
However, with the projected Auger experiment (jCronin 1992 ), the number of detected CR's would 



increase by a factor ~ 50. If £^crit is ~ 2 x 10^*^ eV, as predicted by our model, then a few bright 
sources above lO^^eV should be identified. In addition, for E'er!! = 2 x lO^^eV our model predicts 
~ 10 sources producing more then 5 events at 5 x lO^^eV, compared to ~ 100 such sources 
predicted in a steady source model with a similar Ecrit- 

The observed characteristics of high energy CR sources depend on the bursting rate u and on 
the typical time delay ri only through their product uti. However, and ri may be measured 
separately, if the time delay is either very short, ri < 50yr, or very long, ri ~ lO^yr. In the former 
case, time variability of high energy sources may be detected, while in the latter, which implies 
large magnetic fields, dispersion in CR arrival directions could be measured. The magnetic field of 
our galaxy can also have interesting observable effects: the images of CRB sources should appear 
elongated perpendicular to the direction of the magnetic field, with a predictable correlation of 
the cosmic ray position and energy. For example, a cosmic ray with E = 3 x lO^^eV could be 
deflected by ~ 10° when arriving along the plane of the galaxy. 

The positions of cosmic rays could also be correlated with those of nearby galaxies to see if 
the events producing them occur in normal stellar populations (Waxman, Fisher & Piran T995| ). 



The identification with GRB's could then lead to further constraints on the nature of the objects 
producing these explosions. 
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Fig. 1. — Results of two Monte-Carlo realizations of the bursting sources model with z^ri = 
10~^Mpc~^: Thick solid line- overall spectrum in the realization, shown as the number of cosmic 
rays per unit logE', times the energy. Thin solid line- average spectrum, obtained when the 
emissivity is spatially uniform and not due to discrete sources; notice that this curve is also 
proportional to Dc{E), from eq. (4). Dotted lines- spectra of the five sources having the largest 
CR flux. Short dashed lines- spectra of the five sources that reach the highest fraction of the 
average flux. Filled circles- Fly's Eye data. Open circles- AG AS A data (Icr errors are shown for 
the flux in bins with more than 1 detected events, and for the energy of the highest energy CR's). 
Long dashed lines- the intensity where each experiment should have detected on average one CR 
in each bin, where the bins are equally spaced in log£' (the upper line corresponds to AGASA). 
The normalization of the average flux is chosen to fit the observations at > 2 x 10^^ eV [at lower 
energies, a contribution from iron cosmic rays from other sources is likely to be present (Bird et al. 
1994, Waxman 1995b)]. 
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Fig. 2. — The average number of CR sources as function of energy, for bursting sources model 
with uTi = 10~^Mpc~^ (the dashed hues are added only for visual aid). 



